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Abstract. We include spontaneous symmetry breaking into the functional 
renormalization group equations for the irreducible vertices of Ginzburg-Landau 
theories by augmenting these equations by a flow equation for the order parameter, 
which is determined from the requirement that at each renormalization group (RG) 
step the vertex with one external leg vanishes identically. Using this strategy, we 
propose a simple truncation of the coupled RG flow equations for the vertices in the 
broken symmetry phase of the Ising universality class in D dimensions. Our truncation 
yields the full momentum dependence of the self-energy S(fc) and interpolates between 
lowest order perturbation theory at large momenta fc and the critical scaling regime 
for small fc. Close to the critical point, our method yields the self-energy in the scaling 
form E(fc) = k^a^ {\k\^, |fc|/fcc), where ^ is the order parameter correlation length, kc 
is the Ginzburg scale, and a~{x,y) is a dimensionless two-parameter scaling function 
for the broken symmetry phase which we explicitly calculate within our truncation. 
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1. Introduction 

The functional renormalization group (FRG) has been invented by Wegner and 
Houghton |T] in the early days of the renormalization group (RG) as a mathematically 
exact formulation of the Wilsonian RG. In the past decade this method has gained 
new attention. While there are several equivalent formulations of the FRG involving 
different types of generating functionals, in many cases it is advantageous to formulate 
the FRG in terms of the generating functional F of the one-particle irreducible vertices, 
which can be obtained from the generating functional of the connected Green functions 
via a Legendre transformation [21 [3]. Two different strategies of solving the formally 
exact FRG equation for the functional F have been developed: the first is based on the 
combination of the local potential approximation (LPA) with the derivative expansion 
[H [5]. This approach has been very successful to obtain accurate results for critical 
exponents |H El [6l [7] and is convenient to describe the broken symmetry phase [8], 
because it is based on an expansion in terms of invariant densities, which automatically 
fulfill all symmetry requirements. 

The other strategy, which was pioneered by Morris [3] and has been preferentially 
used in the condensed matter community to study non-relativistic fermions [9l [10], is 
based on the expansion of F in powers of the fields, leading to an infinite hierarchy 
of coupled integro-differential equations for the one-particle irreducible vertices. This 
approach has the advantage of providing information on the momentum- and frequency 
dependence of the vertices. However, there have been only two conceptually different 
attempts to extend the hierarchy of FRG flow equations for the vertices arising from the 
field expansion into the broken symmetry phase. One possibility is to include a small 
symmetry-breaking component into the initial condition for the self-energy fH] [T2| [T3] 
and to check whether this component evolves to a macroscopic value as the RG flow 
is integrated. A disadvantage of this scheme is that the order parameter field and its 
fluctuations do not explicitly appear and that one even has to invest some effort to 
recover simple mean-field results. 

Another possibility to extend the FRG flow equations for the irreducible vertices 
into the broken symmetry phase was proposed in [T3] (see also p[5j). The basic 
idea is to augment the hierarchy of flow equations for the vertices by an additional 
equation for the flowing order parameter, which is obtained from the requirement that, 
at each stage of the RG flow, the vertex with one external leg vanishes identically. 
The purpose of this work is to show how this strategy works in practice. For 
simplicity we shall consider here a simple classical scalar y^^-theory, describing the 
critical behavior of the Ising universality class; generalizations of our method to quantum 
systems are straightforward. For example, this method has recently been used to 
study superconductivity in the attractive electron gas, where the flow equation for the 
order parameter is equivalent to a generalized BCS gap equation including fluctuation 
corrections [T6] . 

Apart from showing how symmetry breaking can be taken into account in the 
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field expansion, we present in this work two additional new results: on the one hand, 
we propose a simple truncation of the exact hierarchy of the flow equations for the 
irreducible vertices in the broken symmetry phase which yields the full momentum 
dependence of the self-energy interpolating between the perturbative regime for 

large momenta k = \k\ and the critical regime for /c — > 0. Different strategies to 
calculate the fc-dependence of the self-energy (and, more generally, the momentum 
dependence of the higher order vertices) has recently been developed in pTj and in 
[T8] . On the other hand, we present in this work an approximate calculation of the 
two-parameter scaling function a~{kC,, k/k^) describing the scaling of the self-energy of 
the system slightly below the critical temperature Tc. Here ^ is the order parameter 
correlation length, and k^ is the Ginzburg scale, which remains finite at the critical 
point p^. The corresponding scaling function {kC, , k / kc) in the symmetric phase 
(i.e., for temperature T > Tc) has recently been discussed in [20]. The fact that the 
Ginzburg scale kc appears in the scaling of thermodynamic variables has been discussed 
in several recent works [21], [22l H], but apparently has been ignored in the older RG 
literature [231 1211 [25]. For models with weak interactions, a universal regime, covering 
the complete crossover from the vicinity of the Gaussian fixed point to the vicinity of 
the Wilson-Fisher-fixed point, exists which can be described completely within a two 
parameter scaling theory [211 [22l H] . For the weakly interacting Bose gas at criticality 
the one-parameter scaling function cr*(fc/fcc) = (y^{oo,k/kc) has been calculated in 
[26| [271 [T71 [18] . However, the full implications of a finite scale kc at criticality and the 
resulting two-parameter scaling theory of the correlation function away from criticality 
has only recently been examined [20] . The field expansion of the FRG allows us to study 
the complete momentum dependence of correlation functions and is thus ideally suited 
to explore the extended universality near the critical temperature T^. 

The structure of the paper is as follows. In Sec. [2] we formulate the exact FRG 
flow equations for the running order parameter Ma and the momentum dependent self- 
energy SA(fc). Guided by the LPA, in Sec.[3]a truncation scheme of the hierarchy of flow 
equations for the irreducible vertices in the broken symmetry phase is introduced, which 
allows us to calculate the momentum dependent self-energy S(fc). In Sec. [Dwe derive 
the two-parameter scaling function for the self-energy using different approximations. 
We end in Sec. [5] with a brief summary and mention further applications of our method. 

2. Exact RG flow equations in the broken symmetry phase 

Our starting point is the following classical action. 



where an ultraviolet cutoff Aq is implicitly understood. In the broken symmetry phase 
the Fourier transform ip^ has a finite vacuum expectation value. 




(1) 



(^fc = + A^k , = (27r)^5(fc)M , 



(2) 
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with (A(/9fc) = 0. By substituting this expression into ([T]) and expanding in powers of 
Aipk, we generate also terms involving one and three powers of Aip^- Demanding that 
the vertex T^l^ associated with the term linear in A(fk should vanish, we obtain the 
magnetization in the Landau approximation, 

" ^ ' -QrUuA, for < " 



This will serve as the initial condition for the flow equation of the order parameter in 
our functional RG approach. 

To derive an exact hierarchy of flow equations for the vertices of our model we 
introduce a momentum cutoff A into the free propagator separating fluctuations with 
small momenta |fc| < A from those with large momenta \k\ > A. Differentiating the 
generating functional F of the one-particle irreducible vertices with respect to A and 
expanding F in powers of the fields, we obtain a formally exact hierarchy of FRG 
flow equations for the vertices 0]. To take into account symmetry breaking, we follow 
the approach proposed in [Hj and demand that for all values of A the flowing vertex 



F^^"* associated with the term linear in the fluctuation A(pk vanishes. This yields a 
renormalization group equation for the flowing order parameter Ma. We would like to 
calculate the true order parameter M = limA^o and the true momentum dependent 
single-particle Green function G{k) = lim\^Q G ^{k) , which we parameterize in terms 
of an irreducible self-energy S(A;), 

For our purpose it is convenient to include the term proportional to in ([I]) into the 
definition of the self-energy. In the broken symmetry phase, the initial condition for the 
self-energy at scale A = Aq is then 

SAo(fc) = rAo + = = -2rAo , (5) 

where Maq is given in ([3]). As we reduce the cutoff, the evolution of the self-energy is 
determined by the following exact RG flow equation [HI [28] . 

1 r d^k' 



5ASA(fc) = -^ I ^GA{k')Tf{k\-k\k,-k) 



^^GA{k')Gj,{k' + k)rf\k, -k - k', k')rf\-k', k + k', -k) 



,(3) 



{dAM^)T)^>{k,-k,0), (6) 



while the flowing order parameter Ma satisfies 

]D 



(9aMa)Sa(0) = - ^/ ^GA(fc)Ff (fc,-fe,0) . (7) 

The flow equations for the three-point vertex F]^'^(fci, k2, k^) and the four point vertex 
F^^(fci, ^2, ^3, ^4) in the presence of symmetry breaking have been written down 
diagrammatically in [H]. In the present work we shall not need these equations. 
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In the broken symmetry phase the proper choice of the cutoff procedure is a dehcate 
matter. The simplest choice is perhaps a sharp cutoff in momentum space, where the 
propagator is for |fc| < Aq given by [3j 

^^^^^ - + ' 

and the corresponding single-scale propagator is 

^'^^^ = -A2 + S^(fc) • 

While in the symmetric phase the sharp cutoff is very convenient [HI [17], it leads to 
technical complications in the broken symmetry phase (see the discussion after fISBl) 
below). These can be avoided if we use a smooth cutoff procedure, which we implement 
via an additive regulator RA{k) in the inverse propagator [5]. The cutoff dependent 
propagator is then 

and corresponding single-scale propagator is 

GA{k) = [-dARA{k)]Gl{k) . (11) 

At this point it is not necessary to completely specify the cutoff function R^^k), except 
that we require it to be of the form [29l |30| 

RAik) = (1 - 4,o)A'^r'^(fcVA') , (12) 

where R{x) is some dimensionless function satisfying R{oo) = and -R(O) = 1. The 
inverse of the flowing wave-function renormalization factor is given by 



dk^ 



(13) 

fc2=0 



where / = — ln(A/Ao). The introduction of Z^^ in ( fT2l) is necessary to preserve the re- 
parametrization invariance of physical quantities under a rescaling of the fields [50] . 
For explicit calculations we shall use the Litim cutoff |31j . 

R{x) = {I - x)e{l - x) . (14) 

Another popular choice is [5] 

— 1 

which has the advantage of being analytic, but leads to more complicated integrals. 



3. Truncated flow equation for the momentum dependent self-energy 

The right-hand side of ([6]) depends on the vertices and T^^^ with three and four 
external legs, which satisfy more complicated flow equations [T^ [28] involving higher 
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order vertices. Keeping in line with the derivative expansion for the effective action [S] 
we truncate the hierarchy as follows 



Ea(0) « ^Ml , (16a) 

rf(fci,fc2,fc3) ~«aMa, (166) 
T^^\ki,k2,k3,ki) ^ ua . (16c) 

This truncation of the field expansion is motivated by the LPA with quartic 
approximation for the effective potential Ues, where one approximates the generating 
functional F by [5] 

T[^] ^ I rf^r[/efrb'(r)], (17) 

with 

UMV]^^[V'-Mlf. (18) 

The completely local character of all correlations in the LPA is known to be a surprisingly 
good approximation in the scaling regime close to criticality [5j. Outside the critical 
regime the LPA fares less well and in general cannot reproduce the structure known 
from perturbation theory (in the case considered here, only the leading order from 
perturbation theory will be reproduced). The LPA combined with a derivative expansion 
converges best if one expands around the local minimum of Ues, see |32]. The condition 
that Ma is the flowing minimum leads to a flow equation for Ma, which in the field 
expansion leads to (I16am6c|) . In this truncation, the exact flow equation ([7]) for the 
order parameter reduces to 

OaMI = -3 1 (0jGA(fc), (19) 

while our flow equation ([2D for the self-energy becomes 

9ASA(fc) = (0jGA(fc') + ^9aM1 - ulMl I ^GA{k')GA{k' + k) 

= - «A / I^^A(fc') - ulMl I ^GA{k')GA{k' + k) . (20) 

By demanding that the flow of Sa(0) is consistent with our truncation ( 116 ap we obtain 
the flow equation for the effective interaction, 

d^UA = -3ul J .^^GA{k)GA{k) . (21) 

The above equations (fT9] - [2T|) form a closed system of coupled integro-differential 
equations for the order parameter Ma, the effective interaction u\, and the momentum 
dependent self-energy SA(fc). In contrast to the LPA, these equations can be used to 
calculate the full fc-dependence of E(fc). Our truncation is similar in spirit but not 
identical to the more elaborate truncation proposed in [18], who also used the LPA as 
a guide to propose a truncation of the hierarchy of flow equations for the momentum 
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dependent vertices generated in the field expansion. However, unlike our equation 
the flow equation for the self-energy proposed by Blaizot et al. |T8] does not involve 
the flowing order parameter, because these authors approach the critical point using an 
expansion around the symmetric state. 

At this point it is convenient to rescale all quantities to reveal their scaling 
dimensions. We define dimensionless momenta q = k/A and the dimensionless coupling 
constants 

Ul =irz5^fA^-V, (22) 

which are considered to be functions of / = — ln(A/Ao). Here Kd is defined by 

(27r)^ 7r^/2r(D/2) ' ^ ' 

where Qp is the surface area of the Z)-dimensional unit sphere. We also define the 
rescaled exact propagator, 

Gi{q) = ^GA{k) 

(25) 



W + rf)(g) + W) ' 

and the corresponding single-scale propagator 

Gi{q)=Ri{q)G]{q), (26) 

where 

= -j-dARA{k) = -(2 - Vi)R{q^) + 2q'i?'(q'). (27) 

Here R'{x) = dR{x)/dx and rji = —di In Zi is the fiowing anomalous dimension. For the 
Litim cutoff R'{x) = —0(1 — x), so that 

Ri{q) = [-2 + r^i{l-q')]eil-q'). (28) 

The rescaled propagator fl25l) depends on the rescaled irreducible self-energy, 

(q) = §s^(fc) • (29) 

By construction, the constant part of the rescaled self-energy is 

rf)(0) = |Mf = ||Ml. (30) 

The flow of rp''(0) is thus determined by the flow of and m/, which in our truncation 
is given by 

diMf = {D-2 + r^i)Mf + 3 I Gi{q) , (31) 

Jq 

diui = {A-D- 2rii)ui + ?>ui [ Gi{q)Gi{q) , (32) 

J a 



(36) 
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where = Vf^ J d^q. To calculate 77^, we need the momentum dependent part of the 
rescaled self-energy, 

li{q)=T^P{q)-Tf\Q), (33) 

which satisfies 

dai{q) = (2-rii-q- V^Mq) + 7^(9) , (34) 

where 

7Kq) = u^M^ I Gi{q')[Gi{q' + q) - G,{q')] . (35) 

Jq' 

The flowing anomalous dimension is then given by 

^'=^ 

The reason why using a sharp cutoff in the broken symmetry phase leads to technical 
complications is that in this case the expansion of the integral ( l35l) for small q would 
start with a non-analytic term proportional to \q\, which requires the introduction of 
an additional relevant coupling constant [17] . On the other hand, with the Litim cutoff 
f ll4p or the analytic cutoff (1151) the leading term in the expansion of '^i{q) is proportional 
to q^. In this case the expansion of the right-hand side of fl35l) for small q yields for the 
flowing anomalous dimension 

m = -uM j G,{q)[G]{q)[l + R'{q')] + ^[2G}{q)B!' {q') - 4Gf(q)[l + F^{q')f]] , 

(37) 

where E!\x) = <fR{x)/dx^. For the Litim cutoff E!'{x) = S{1 — x). 

4. Two parameter scaling in the broken symmetry phase 
4.I. Truncation with only marginal and relevant couplings 

In the simplest self-consistent approximation, we expand the momentum dependent part 
7;(q) of the two-point vertex on the right-hand side of our flow equations (l3Tl [32l E 
to first order in q^. Since by definition 

dVPiq 



1 - 



^ _ d-fi{q) 



dq 

q2=0 



(38) 

q2=0 



this amounts to approximating the propagator on the right-hand side of the flow 
equations (EH [321 [M]) by 



where 



Pi = T?\0) = ^Mf, (40) 
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Figure 1. Typical non-critical flow of the coupling parameters and ui obtained 
from (ETl) and m in £> = 3. 



see ( 130!) . The resulting system of flow equations for the couphng constants and 
ui together the flow equation diZi = —rjiZi for Zi are equivalent to the quartic 
approximation for the effective potential with wave-function renormalization [5]. Using 
the Litim cutoff dHD, equations ( l3Ti) and ( l32i) become 

ml = {D-2 + ^,)Mf - ^%~lj^~^f GKO), (41) 

d,u, = {A-D- 2r,,)u, - ^%\^^^^^ <Gm. (42) 

where G;(0) ~ [1 + p/]~^ is the rescaled propagator at zero momentum. Moreover, with 
the Litim cutoff the flowing anomalous dimension ( 1371) is simply 

VI = ^uMG'M. (43) 

Equations pTT - HSl) form a closed system of differential equations for Mf, ui and rji which 
can easily be solved numerically. To find the flow along the critical surface, we need to 
fine tune carefully the initial values mq and Mq . A typical flow of the rescaled parameters 
as a function of / is shown in figure (H while in figure [2] we show the flow schematically in 
the (u/, M;^)-plane. The Wilson-Fisher fixed point in D = 3 is in this approximation at 
~ 0.942 and Mf ^ 1.022. As can be seen in figure [2l the couplings initially flow very 
slowly and stay close to their initial values in the vicinity of the Gaussian fixed point. 
At a characteristic scale Ic they are rapidly attracted by the Wilson-Fisher fixed point 
where the flow is again almost stationary. Finally, at the scale all non-critical RG 
trajectories rapidly move away from the Wilson-Fisher point and the /-dependence of the 
couplings Ui and Mf is determined by their scaling dimensions, ui oc e^', Mf oc e^^"^-*', 
where e = 4 — the flow of the un- rescaled couplings ma and M\ then stops. 

What determines the two characteristic scales Ic and The momentum scale 
kc = Aqc"''' associated with Ic measures the size of the Ginzburg critical region. For 
small initial values uq the logarithmic scale Ic is given by [221 [El IlH] 

lc^-ln(^), (44) 



e \uo 
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Figure 2. Qualitative flow diagram for the couplings ui and M^. The arrows sketch 
a nearly critical flow, with the size of the arrows representing the velocity of the flow. 
The dots mark the Gaussian flxed point (G) and the Wilson-Fisher fixed point (WF). 



where is the value of ui at the Wilson-Fisher fixed point. This scale can be derived 
from f HTl) and ( H2i) by approximating Gi{0) ~ 1 and t]i ^ 0. In the intermediate regime 
h '^l the fiowing can then be replaced by a constant ^ M'^ = 6/[D{D — 2)], 
while for all I the solution of (H2l) can be approximated by 

^ ^ I (45) 

where = De/6. The numerically obtained flow shown in flgure [1] further reveals that 
the scale Ic is also characteristic for the /-dependence of M^. 

Non-critical flows which describe the system at T < Tc eventually obey Mf /M^ 3> 
Gf{Q) and u^jui ^ Cf (0). In that case, the solution for Mf and ui depend exponentially 
on /, corresponding to trivial scaling. The unrescaled variables u\ and M\ then approach 
flnite limits and also the physical correlation length ^, which is deflned via 

r' = lim[ZASA(0)], (46) 

remains flnite. The scale in flgure [1] is related to ^ via = Aqc"'* or equivalently 
2Z, = -ln[lime-2'rf)(0)]. (47) 

Z^oo 

From the linearized flow around the Wilson-Fisher flxed point we obtain in D = 3 
the critical exponents u ~ 0.553 and t] ^ 0.099. Very similar results are obtained using 
the analytic cutoff (ITSl) . The poor comparison of our results to the established values 
u = 0.64 and rj = 0.044 (see [5]) can be traced to the low order truncation of our 
effective potential f|T8|) . see [32l [33]. Close to D = 4 we obtain in the leading order 
e-expansion r] ~ e^/12 and z/~l/2-|-e/12. While the result for u is correct, the value 
for Tj compares badly with the known expansion t] ~ e^/54. This clearly shows the 
limitations of a low-order effective potential approximation. 

4-2. FRG enhanced perturbation theory 



So far, we have truncated the self-energy retaining only its marginal and relevant parts. 
This is a good approximation for small momenta k. On the other hand, for large k this 
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approximation cannot correctly reproduce the momentum dependence of the self-energy 
which arises from perturbation theory. To leading order in the relevant dimensionless 
bare coupling uq = mao^^~^; the perturbative momentum dependence of in the 
broken symmetry phase is given by [20] 



with 



where 



S(fc)-S(0)=r'Aoro-(A;0, (48) 
Aa,{x) = ^-^[x{0)-x{x)]+Oiul), (49) 

^^^^ = / (271)^ [p'^ + l][(p' + py + ly ^^0) 

( l48] - [50l) are only accurate sufficiently far away from the critical point where ^ and the 
relevant dimensionless coupling uq = maq^'^"^ are small. 

We now present an improved approximation for the momentum dependent self- 
energy which we call FRG enhanced perturbation theory since it embeds the perturbative 
expansion into a functional renormalization [31] such that it reproduces exactly the 
leading order perturbative behavior for large k. However, in contrast to perturbation 
theory, FRG enhanced perturbation theory does not suffer from any divergencies; it 
yields an explicit description of the entire crossover to the critical regime and gives 
reasonable results even at the critical point. 

Quite generally, the physical self-energy can be written as an integral over the entire 
RG trajectory [20] . 

S(fc) - S(0) = / dle-^'+^o'i-rir^le'k/Ao), (51) 
Jo 

where S(0) = Z^^^~'^. In general, the expression for the inhomogeneity 7(g) will 
depend also on the momentum dependence of the three- and four-point irreducible 
vertices, as can be inferred from ([6]). To make progress, we employ the truncation 
f ll6atll6c|) which leads to the approximation ( l35l) for the inhomogeneity 'y{q)- Only the 
momentum independent parts of the three- and four point vertices enter and 7(9) is 
then completely determined by the self-energy and the order parameter alone. While 
this greatly simplifies the calculation of the self-energy since it leads to a closed set of 
equations, solving ( l5Ti) remains non-trivial since the solution for the self-energy must be 
determined self-consistently. The FRG enhanced perturbation theory provides for a non- 
self-consistent approximation to the solution of (15T!) . In the FRG enhanced perturbation 
theory, the calculation of the subtracted inhomogeneity 7(9) via (1351) is simplified by 
keeping only the first two terms in a momentum expansion of the self-energy. This 
amounts to the substitution 

rp)(9)^rf)(0) + (l-Z,)q^ (52) 

and the approximation (139!) for the propagator on the right-hand side of (135|) . Within 
this approximation, the fiow of M^, ui, and rji are determined from (HTI - 1431) . A similar 
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truncation strategy has been adopted in [T7] for the symmetric phase of the 0(2)-model, 
and in [35] to calculate the spectral function of the Tomonaga-Luttinger model. A 
comparison with the completely self-consistently determined self-energy is presented at 
the end of this section where we show that the error arising from the non-selfconsistency 
of the FRG solution is extremely small. Perturbation theory is recovered when the flow 
of the running couplings is approximated by their trivial /-dependence arising from their 
scaling dimensions. The self-energy can now be expressed in terms of a two-parameter 
scaling function, 

^{k)=kla-{x,y), (53) 

with X = kC, and y = k/kc- The ratio of these variables is then x/y = e'"'"'"'. If we 
introduce 

Aa-{x,y) = a-{x,y) - kj^O) = a' {x,y) ~ y^ / Zx\ (54) 

this leads to 

poo ^ poo 

Aa-ix,y) = / dZe-2('-'^)+/o'^-^^^,(e'-'c^) = / dpp-'Z-l,^^^/^^%^,^(^,/y){p), (55) 

Jo Jye-'c 

where we substituted p = ye^'^" and used Zi = e~ -^o '^'^^^ . The asymptotic behavior of 
7j(q) for small q is 7i(q) ~ ViQ'^- large q it approaches a constant which, using the 
Litim cutoff, is 

Jim 7^9) - 2«'^'Gi^(0)%^^^ • (56) 

In D = 3 the function 'yi{q) can be calculated analytically for the Litim cutoff; the 
result is given in the Appendix. At the critical point x ^ oo since the correlation 
length diverges, so that the scaling function reduces to 

Aa' {oo, y) = a' {oo, y) = a^{y). (57) 

The asymptotic behavior of the scaling function for both very small and very large y 
follows directly from fl55|) . For y <^ 1, i.e. in the critical long wavelength regime |20j . 
the lower limit of integration may be replaced by zero and all coupling parameters may 
be replaced by their fixed point values. Then we find 

a^y) ^ Any'-", (58) 

where rj is the fixed point value of rji and 

POO 

Ad= / dpp^-^M. (59) 



with 7*(p) = lim/^oo 7; (p) • In D = 3 we obtain numerically ^3 1.075. In the critical 
long wavelength regime y ^ 1 one may approximate all couplings by their initial values. 
The scaling function then approaches the constant value 

aM ^ ^ulMlCm. (60) 

Such a constant plateau is expected from the structure of the truncation employed. 
In fact, for k ^ kc one expects that the momentum dependence of the self-energy is 
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In y=ln (k/k^,) 



Figure 3. Typical beliavior of the two-parameter scaling function Act {x, y) defined 
in (|55p for x = e^'^^^y. The initial coupling parameters are uq — 0.005 and 
= 1.9482092, which yields Ic 5.23 and U ^ 11.29. 

that of lowest order perturbation theory, see the discussion at the beginning of this 
subsection. However, an effective ultra-violet cutoff is now provided by which 
regularizes the theory in place of the correlation length which is infinite at criticality. 
While this is indeed the leading order correction to the self-energy in an expansion 
in powers of the bare interaction strength, the correct form of the self-energy should 
further display a \n{k/kc) dependence at large k with a pre-factor which is quadratic in 
the bare interaction [TTl [20] . The reason for the absence of such a term in the present 
approximation is that our truncation for the four-point vertex in fll6cp does not take 
vertex corrections into account. In the macroscopically ordered regime fc.^ <^ 1 we find 
that — S(0) vanishes as A;^, as can be seen in figure [3l 

For the present model, it is of course possible to calculate the solution of the coupled 
integro-differential equations (fT9] - [2T!) exactly without any further approximations. One 
obtains completely self-consistent solutions for the flow of the self-energy SA(fc) and the 
order parameter if the truncation of the momentum dependence of Sa and G\ on the 
right-hand sides of (|T9] - [2Ti) using the substitution ( l52l) is omitted. A comparison of the 
completely self-consistent solution a^^{(X),y) for the scaling function with the scaling 
function cr^(oo,?/) obtained within the FRG enhanced perturbation theory, i. e. with 
the help of the substitution (1521) on the right-hand sides of (fT9] - [2Tl) . is shown in figure 
m Obviously, the relative error due to the substitution fl52|) is remarkably small, so that 
we conclude that our substitution ( l52i) is quite accurate. 

Of course, one could easily improve the approximations presented here. A 
straightforward extension would be to truncate the effective potential Ue?i in (fTTll at 
some higher order n > 2, 

n (m) 
m=2 ^ 

To arrive at the flow equations of the parameters . . . , one would need to 
take into account the flow equation of the lowest n vertices of a field expansion. We 
have done so up to n = 5; a significant improvement of the critical exponents rj 
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Figure 4. Comparison of two different approximations of the self-energy at criticality, 
as discussed in the text: self-consistent numerical solution a^^{oo, k/kc) of p9H2T|) and 
solution a~{oo,k/kc) based on the substitution ((52| on the right-hand sides of these 
flow equations. 

and u in the Ising model would however require approximately n = 10, as is known 
from previous investigations of the derivative expansion [33]. Note that our FRG 
enhanced perturbation theory can also be adopted if arbitrarily higher orders in n are 
included. This is because within our approach the flow of the coupling constants which 
parameterize the local potential depends only on the lowest order momentum expansion 
of the self-energy and can be calculated exactly as in the usual derivative expansion |33j . 
Once the flow of the local potential is known, the higher order momentum dependence 
of the self-energy can be determined. Of course, not all models have a structure as 
simple as the one discussed here which permits to include arbitrarily high orders of the 
local potential. One may wonder whether our approach is also useful to describe more 
complicated systems. It is certainly expected to be useful in non-critical interacting 
systems, as e. g. interacting bosons in two or three dimensions, where a low order 
truncation of the effective action should suffice [36]. An accurate description of the 
momentum dependent self-energy of more complicated and possibly critical systems, 
such as e.g. frustrated spin models [6], is a very challenging task which we have not yet 
attempted. While including all orders of the local potentials would then be prohibitive, 
a low order truncation might yet be qualitatively correct. 

5. Summary and conclusions 

Let us summarize the three main results of this work: 

(i) We have demonstrated how symmetry breaking can be included in the exact 
hierarchy of FRG flow equations for the irreducible vertices within the framework 
of the field expansion, using the approach from [H]. The basic idea is to require 
that the vertex with one external leg vanishes identically for all values of the fiow 
parameter, which yields an additional fiow equation for the order parameter. Our 
method differs from that employed in [111 [I2l [13], which do not explicitly include 
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a flow equation for tlie order parameter. 

(ii) Guided by tlie LPA, we liave proposed a simple truncation of the FRG flow 
equation for the momentum dependent self-energy which yields a reasonable 
interpolation between the perturbative regime for large momenta and the critical 
regime for k —>■ 0. We have also pointed out that a sharp 0-function cutoff (which 
is still very popular in the condensed matter community) is not suitable to analyze 
the broken symmetry phase, because it generates a non-analytic term proportional 
to |fc| in the two-point vertex for momentum scales below the RG cutoff A 

(iii) Using the above truncation, we have calculated the two-parameter scaling function 
a^{kC,, k/kc) describing the scaling of the self-energy in the broken symmetry phase 
slightly below the critical temperature. Similar to the case T > T^. discussed in [20] , 
the scaling function depends on two parameters, involving the correlation length 
1^ and the Ginzburg scale kc. The latter remains finite at the critical point and 
measures the size of the critical region. 

It is straightforward to generalize the method described here to study spontaneous 
symmetry breaking in quantum mechanical many-body systems. For example, our 
approach can be used to obtain fluctuation corrections to the BCS gap equation in the 
attractive Fermi gas [16]. Note that beyond the BCS approximation the superconducting 
order parameter should be distinguished from the off-diagonal self-energy associated 
with the single-particle Green function, so that it is important to introduce both 
quantities into the FRG as independent parameters. We are currently using our method 
to study the interacting Rose gas [36], where a truncation similar to the one discussed 
here yields corrections to the Bogoliubov mean-field approximation for the diagonal and 
off-diagonal self-energies which are consistent with the Hugenholtz-Pines theorem |37j . 
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Appendix: Analytical form of 7/(q) in three dimensions 

Using the Litim cutoff [31] given in ( |T^ we can rewrite the function 7z(q) defined in 
([35D as 




(A.l) 
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where the functions 7^^^ (g) and •y^' (g) are defined by 

Pi + il + q?' 



M^ufGf{0)\A{q)+Bi{q) In 
-Ci{q) 



arctan 



Pi + l 



— arctan 



q 



(A.3) 



with 

Ai{q) : 

Biiq) : 
Ci{q) 

Diiq) 

EM) 



+Q{q) 



A + (1 + q? 



arctan 



Pi + (1 - 
+ arctan 



q)' 
l + q 



(A.4) 



^|60[1 + pi][A - r]i{l + pi)] - 30g[4 - r]i + Api{3 - 2r]i) - Ir^ipJ 



16g 



{4[g2 - - 1] + r^^pf ^ 2p,(l - 3^^) + (g^ _ i)2]|^ 



I -^[..(5 + P.) 

1 + P; 
16g 



:i+pO[i-|(2+3p, 



jr]/ - 4 + pi[r]i{2 + p^) - 4] - 2g2[3r/zpz + r/z - 2] + r/^g^j. 



(A.5) 
(A.6) 

(A.7) 
(A.8) 
(A.9) 
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